Role of Boron p-Electrons and Holes in Superconducting MgB 2 , and other Diborides: 
A Fully-relaxed, Full-Potential Electronic Structure Study 
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We present the results of fully-relaxed, full-potential electronic structure calculations for the new 
superconductor MgBi, and BeBi, NaB2 and AIB2, using density-functional-based methods. Our 
results, described in terms of (i) density of states (DOS), (ii) band-structure, and (iii) the DOS 
and the charge density around the Fermi energy Ef, clearly show the importance of B p-band for 
superconductivity. In particular, we show that around Ef, the charge density in MgBi, BeBi and 
NaB2 is planar and is associated with the B plane. For BeB^ and NaB2, our results indicate 
qualitative similarities but significant quantitative differences in their electronic structure due to 
different lattice constants a and c. 
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The intricate interplay of interactions between elec- 
trons, holes and phonons is quite evident in the recently 
discovered superconductor Magnesium Boride, MgB 2 §. 
The superconducting temperature of around ~ 39 K in 
MgB 2 is almost at the extreme end of T c 's predicted by 
the BCS theory for conventional superconductors. The 
ongoing experimental investigations suggest it to be a 
BCS type superconductor with confirmed isotope effects 
Pj and other related effects ||. However, more exper- 
iments are needed to completely characterize the na- 
ture of interaction responsible for superconductivity (SC) 
in MgB 2 . Theoretical studies of MgB 2 , since the an- 
nouncement of SC by Akimitsu, have concentrated on (i) 
the nature of interactions responsible for SC using model 
Hamiltonians and (ii) the electronic structure |^-||] 
and the linear response |,|||§ of normal state MgB 2 
using ab initio methods but with the experimentally ob- 
tained lattice constants a and c. 

In order to understand SC in MgB 2 , be it due to heav- 
ily dressed holes Q , hole-doped covalent bonds §| or sim- 
ply due to electron-phonon interaction |ic[ | , it is essential 
to have an accurate electronic structure description of 
the normal state MgB 2 . It is also becoming clear that 
the boron layer and the two-dimensional cr band due to 
Px.y orbitals in that layer hold the key to SC in MgB 2 
(we present evidence of that) fj-g]. Thus an understand- 
ing of how the a band and the p x ,y orbitals respond to 
changes in the chemical environment, lattice constants 
etc. could be quite useful in understanding SC in MgB 2 
and in synthesizing new materials with MgB 2 -\ike SC. 
Fortunately, these changes around B layer can be mim- 
icked without destroying its two-dimensional character 
which is essential for SC. For example, the effects of re- 
moving one electron and adding one electron around the 
B layer can be studied by considering NaB 2 and AIB 2 , 
respectively. Similarly, a study of BeB 2 (and NaB 2 ) will 
provide the effects of changing the lattice parameters c 
and to a lesser extent a. Thus an accurate and reliable 



systematic study of the electronic structure of MgB 2 , 
BeB 2 , NaB 2 and AIB 2 , coupled with the experimental 
facts of SC in MgB 2 and no SC in AIB 2 , can provide sub- 
stantive clues about the nature of interaction responsible 
for SC as well as help in synthesizing new superconduct- 
ing materials with T c 's comparable to that of MgB 2 . 

As the nature of bonding in boron layer depends crit- 
ically on lattice constants a and c, and since for SC one 
needs the electronic structure within a small energy inter- 
val of ±hu>D ( = Debye frequency) around the Fermi 
energy, which gets affected the most due to change in 
lattice constants, it is of utmost importance to theoreti- 
cally optimize the lattice constants a and c. Some of the 
earlier studies of MgB 2 , which are based on the experi- 
mental lattice constants, have provided a good overall de- 
scription of the electronic structure | ||jg , |ll| , |l2[ and some 
others have used the Fermi surface [|6| to discuss SC in 
MgB 2 . Previous study jjj of BeB 2 was carried out with 
assumed values for the lattice constants a and c. A re- 
cent calculation || of BeB 2 gives optimized values for a 
and c but we find these values to be not fully converged. 
In our opinion, there is a clear need to include (a) struc- 
tural relaxation of the lattice, (b) a detailed description 
of electronic structure in a small interval around Fermi 
energy, and (c) a comparison with other related com- 
pounds which have either the potential of becoming an 
MgB 2 -\ike superconductor or help in understanding SC 
in MgB 2 , in a density-functional-based electronic struc- 
ture calculations of MgB 2 . The present work is a step in 
that direction. 

In the following, we present the results of fully-relaxed, 
full-potential electronic structure calculations for MgB 2l 
BeB 2 , NaB 2 and AIB 2 in PQ/mmm crystal structure, 
using density-functional-based methods. We analyze our 
results in terms of (i) density of states (DOS), (ii) band- 
structure along symmetry directions, and (iii) the DOS 
and the electronic charge density in a small energy win- 
dow around the Fermi energy, Ep, in these systems. 
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For our calculations, we have used ABINIT code [ p^[ , 
based on psuedopotentials and planewaves, to optimize 
the lattice constants a and c of MgBi 1 BeB2, NaB2 and 
AlBi . In order to study the band-structure along sym- 
metry directions and site- and symmetry-decomposed 
densities of states, we have calculated, self-consistcntly, 
the electronic structure of these compounds using our 
own full-potential program as well as LMTART package 
with the optimized lattice constants. For studying 
the charge density in a small energy window around Ep 
we have used the Stuttgart TB LMTO package [HJ . 

Based on our calculations, described below, we find 
that in MgB-i (i) B p-band is crucial for SC due to its 
proximity to Ep. In particular, it is evident that within 
a small energy window around Ep, the charge density is 
two-dimensional and it is associated with the B plane, (ii) 
structural relaxation leads to significant changes around 
Ep in the densities of states and charge density. For the 
hypothetical compounds BeB2 and JV0B2, our results in- 
dicate qualitative similarities but significant quantitative 
differences in electronic structure due to very different 
lattice constants a and c of these compounds. A com- 
parison of band-structure of BeB2, NaB^, MgB^ and 
AIB2 along symmetry directions reveals that the point 
moves down, starting from above Ep in NC1B2 to well be- 
low Ep in AIB2- If SC is due to holes then one would 
expect BeB2 and NaB2, if they can be synthesized, to 
become superconducting but not AIB2- 

Before describing our results in detail, we pro- 
vide some of the computational details of our cal- 
culations. The structural relaxation was carried out 
by the molecular dynamics program ABINIT with 
Broyden-Fletcher-Goldfarb-Shanno minimization tech- 
nique (ll| using Troullier-Martins pseudopotentials jlfj] 
, 512 Monkhorst-Pack jn| k-points and Teter param- 
eterization Jl3| for the exchange-correlation. The ki- 
netic energy cutoff for the plane waves was 110 Ry. The 
charge self-consistent full-potential LMTO calculations 
were carried out using 2K-panels with the generalized gra- 
dient approximation for exchange-correlation as given by 
Perdew et al. Ilq] and 484 k-points in the irreducible 
wedge of the Brillouin zone. The basis set used consisted 
of s, p, d and / orbitals at the Mg site and s, p and 
d orbitals at the B site. The potential and the wave- 
functions were expanded up to l rnax = 6. The input to 
the tight-binding LMTO calculations, carried to charge 
self-consistency, were similar to that of the full-potential 
calculations except for using spherically symmetric po- 
tential and the space-filling atomic spheres. 

In Table I we compare the calculated lattice constants 
a and c for MgB2, BeB2, NaB2 and AIB2 with the ex- 
perimentally observed values [ ^9|j20| , indicated in the 
parentheses. The calculated lattice constants for MgB2 
and AIB2 are smaller than their respective experimen- 
tal values by roughly 1%. It compares well with the 
differences mentioned in 101 for the lattice constants of 



MgB<2- For BeB2 and NaB2 the values obtained for c/a 
are 0.98 and 1.43, respectively. Although, one expected 
the change in c/a ratio to be small in £?e£>2 and large 
in NaB2 as compared to MgB2, however, the c/a value 
for NaB2 is larger than expected. Because of the large 
difference in the c/a ratio of BeB2 and NaB2 they can 
be used to explore the importance of inter-layer coupling 
for SC in these systems. A recent calculation Q of BeB2 
gives a somewhat smaller value for a and c. We have used 
the lattice constants as given in Ref. j|] to calculate the 
total energy of BeB 2 using the full-potential LMTART 
program, and we find that the lattice constants as given 
in Table I of our paper correspond to lower energy by 
about 3 mRy. We have also run the ABINIT code with 
the information given (which is not complete) in Ref. ||] 
to optimize the lattice constants of BeB2- The results 
are in agreement with the values given in Table I. In Ta- 
ble I, we also show the plasma frequencies and the total 
density of states per spin, n(Ep), at Ep for these com- 
pounds. The plasma frequencies along x {y ) direction 
are very similar for all four diborides but along the z di- 
rection it correlates with the c/a ratio of the compound. 
A smaller c/a ratio leads to enhanced coupling in the z 
direction. 

We show in Fig. 1 the total density of states for BeB2, 
NaB 2 , MgB 2 and AIB 2 calculated using the ABINIT 
program at the optimized lattice constants of these com- 
pounds as given in Table I. The gross features of the DOS 
of the four compounds are similar if one takes into ac- 
count the differences in the total number of valence elec- 
trons. The bottom of the band is the deepest for AIB 2 
which has a total of nine valence electrons and shallow- 
est for NaB 2 with only seven valence electrons. However, 
the isoelectronic structures BeB2 and MgB2 show sub- 
stantial differences in their DOS near the bottom of the 
band. This is due to the fact that the smaller c/a ra- 
tio leads to enhanced repulsion which pushes the s and 
p electrons at the Mg site downward and, at the same 
time, diminishes the DOS for B s and p electrons in the 
middle of the band. In the case of AIB2 we find that the 
B p-band is completely inside Ep. 

An important factor in determining the superconduct- 
ing temperature in conventional superconductors is the 
DOS within an interval of ±Iiwd at Ep, In order to an- 
alyze our results in detail we have plotted in Fig. 2 the 
symmetry-decomposed DOS at the B site within a small 
energy interval around Ep. From the DOS for AIB2, it is 
evident that it will not be superconducting, not at least 
in the same sense as that of MgB 2 , because the bands 
with x and y symmetry are completely filled. According 
to our calculations as shown in Fig. 2, NaB2 is likely to 
show SC with enhanced T c . 

To see how the differences in the lattice constants and 
the number of valence electrons affect the band-structure, 
in Fig. 3 we have plotted the band-structure along the 
symmetry directions for BeB 2 , NaB2, MgB% and AIB2- 
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We find that these differences lead to significant changes 
in the band-structures of BeB2, NaB2, MgB^ and AlB^- 
A measure of these changes is evident in the movement 
of r$ point with respect to Ep as one goes from BeB-2 to 
AIB2. Due to different number of valence electrons, the 
T§ point is well above Ep in NaB2 but well below Ep 
in AIB2- The differences between the band-structures of 
BeB2 and MgB2 arise due to differences in the lattice 
constants. Our calculations show that r$ point is inside 
Ep by about 5 mRy in BeB2, while it is slightly above 
Ep in Mg£>2- It is interesting to note that the location 
of T$ point, around which the Fermi surface in MgB^ 
is hole- like f|,|2l|], has been correlated with SC in MgB 2 
and other diborides Our Z-character analysis shows 
substantial p x and p z characters along T-A direction near 
Ep in all the diborides. 

To complete our analysis we have also calculated the 
electronic charge density within a 5 mRy energy win- 
dow around Ep for BeB2, NaB2, MgB2 and AIB2 com- 
pounds using the TB LMTO method with the lattice con- 
stants as given in Table I. In these calculations the atomic 
sphere radii were adjusted to minimize the discontinuity 
in the Hartree potential across the atomic spheres which 
have been found to give reliable results ||^3| . To reduce 
the overlap errors in NaB2, which has a large c/a ratio, 
we have inserted an empty sphere. In Fig. 4, we show the 
charge density calculated within a 5 mRy energy window 
around Ep for these compounds. For M3B2, BeB2 and 
NaB2we find very similar charge density, which is due to 
p x>y orbitals in the B plane. It should be noted that the 
holes in MgB2 also originate from these planar orbitals in 
the B plane. In contrast, the charge density around Ep 
in AIB2 has a three-dimensional character and it consists 
mostly of Mg s- and B p-electrons. 

In conclusion, we have presented the results of our 
fully-relaxed, full-potential electronic structure calcula- 
tions for BeB2, NaB2, MgB2 and AIB2 in PQ/mmm 
crystal structure using density-functional-based meth- 
ods. We have analyzed our results in terms of the density 
of states, band-structure, and the DOS and the electronic 
charge density around Ep in these systems. For MgB2, 
we find that the p-band of B is crucial for SC due to 
its proximity to Ep. In particular, we have shown that 
within a small energy window around Ep, the charge den- 
sity is two-dimensional and it is associated with the p x>y 
orbitals in the B plane. For BeB2 and NaB2, our results 
indicate qualitative similarities but significant quantita- 
tive differences in the electronic structure due to differ- 
ences in the lattice constants. For example, we find that 
the r§ point in BeB2 is just below the Fermi energy while 
in NaB2 it is well above the Fermi energy. 
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TABLE I. The calculated lattice constants a and c, the 
plasmon frequencies, and the density of states at the Fermi 
energy. The experimental lattice constants are shown in the 
parentheses. 





a (a.u.) 


c (a.u.) 


(eV) 


<4 ( eV ) 


n(E F ) 


BeB 2 


5.49(5.558) 


5.41(5.425) 


8.1 


11.8 


3.29 


NaB 2 


5.70 


8.15 


7.75 


3.40 


6.75 


MgB 2 


5.76 (5.834) 


6.59 (6.657) 


7.04 


6.77 


4.70 


AIB 2 


5.63 (5.681) 


6.13 (6.149) 


8.11 


10.02 


2.66 
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FIG. 1. The total density of states calculated at the opti- 
mized lattice constants using the ABINIT program. 
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FIG. 3. The band-structure along symmetry directions 
calculated at the optimized lattice constants using the 
full-potential LMTO method. 
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FIG. 2. The partial B density of states around Fermi en- 
ergy calculated at the optimized lattice constants using the 
full-potential LMTO method. 



FIG. 4. The isosurfaces of charge density within a 5 mRy 
energy window around the Fermi energy in the primitive cell 
of BeB 2 (top left), NaB 2 (top right), MgB 2 (bottom left) 
and AIB 2 (bottom right), calculated at the optimized lat- 
tice constants using the TB LMTO method. The values of 
the isosurfaces (8,6,4,2,1) decrease as one moves away from 
the B plane. The unit for charge density in the volume is 
10~ 4 [e/(a.u.) 3 ] and the view is off-diagonal at an angle of 
30°. 



